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ABSTRACT 

The Crab Nebula has beeome established as the standard eandle for TeV gamma-ray astronomy using 
the atmospherie Cherenkov teehnique. No evidenee for variability has been seen. The spectrum 
of gamma rays from the Crab Nebula has been measured in the energy range 500 GeV to 8 TeV 
at the Whipple Observatory by the atmospheric Cherenkov imaging technique. Two methods of 
analysis involving independent Monte Carlo simulations and two databases of observations (1988-89 
and 1995-96) were used and gave close agreement. Using the complete spectrum of the Crab Nebula, 
the spectrum of relativistic electrons is deduced and the spectrum of the resulting inverse Compton 
gamma-ray emission is in good agreement with the measured spectrum if the ambient magnetic field 
is about 25-30 nT. 

1. INTRODUCTION 

The Crab Nebula has become a standard candle in TeV astronomy; it has been detected by many 
groups and its integral flux appears constant. (We refer to the “steady” emission, not emission modu¬ 
lated at the pulsar period. We have not detected the latter at TeV energies. See paper by G. Gillanders 
et ah, in these proceedings.) It is also well on the way to becoming a standard candle with regard to 
TeV spectral content. Synchrotron emission from the Crab covers a remarkably broad range termi¬ 
nating at about 10* eV, where a new component attributed to inverse Compton scattering begins. It is 
this component that we detect. The TeV spectrum is sensitive to the primary electron spectrum, the 
nebular magnetic field and the spatial distribution of electrons and magnetic field within the nebula. 

In this paper we briefly describe two methods for extracting TeV spectra, compare results from 
the Whipple Observatory Imaging Cherenkov Telescope for the 1988/89 and 1995/96 observing sea¬ 
sons and comment on implications for the physics of the nebula. The methods, developed at Iowa 
State University and the University of Eeeds, are based on Monte Carlo simulations using completely 
independent code and use different approaches in determining the overall gain of the Cherenkov tele¬ 
scope. In the ISU approach, the gains of the photomultiplier tubes, mirror reflectivities, etc., are 
measured and combined to find the overall gain. In the Eeeds approach, the observed brightness 
distribution of cosmic-ray images is combined with cosmic-ray simulation results to determine the 
overall gain of the telescope. The methods are described in detail in “Paper I,” Mohanty et ah, (1997), 
and the resulting TeV spectrum is put into context of other observations with implications for the 
physics of the nebula in “Paper II,” Hi 11 as et ah, (1997). The latter paper also compares our Crab TeV 
spectrum with those of other groups. 


A straightforward approach to the determination of TeV spectra was developed at Iowa State Uni¬ 
versity. Three components are required. First, a method of distinguishing gamma-ray images from 
background cosmic-ray images. The standard method is to use “supercuts” as described in, e.g.. 
Punch et al. (1996). The images are characterized by second order moments giving the width, length, 
distance of the image centroid from the optic axis and alpha, the angle by which the image major axis 
misses passing through the optic axis. More that 99% of the background can be rejected by requiring 
small values of width, length and alpha. However, this procedure results in a strong bias against the 
images of higher energy gamma rays which tend to be longer and broader and hence more cosmic-ray 
like. Mohanty (1995) has modified the procedure so that the image selection criteria depend on the 
total brightness or size of the image as well. (The size can be used as an estimate of the energy.) 
The telescope collection areas for simulated gamma rays for the 1995/96 season using standard and 
“extended” supercuts is shown in Fig. 1. 

The second component needed is 
a way to estimate the energy of each 
gamma-ray image. Two desirable crite¬ 
ria are (a) good resolution and (b) negli¬ 
gible bias. The former is important to de¬ 
tect small structures in the spectrum and 
the latter is important to avoid distortions. 

We obtained a resolution of tsE jE ~ 0.36 
with negligible bias by using a second 
order polynomial in size and distance as 
described in Paper I. The energy resolu¬ 
tion function is, to a good approximation, 

Gaussian in the variable log(£’). It is plot¬ 
ted at an arbitrary energy in Fig. 1. 

In data taken for spectral analysis, 
each on-source observation is followed by 
an “off-source” observation covering the 
same range of range of azimuth and ele¬ 
vation angles. The images for both types 
of observation are selected for gamma-like events. The estimated energies from corresponding obser¬ 
vations are histogrammed and the difference ascribed to gamma rays from the source. This can then 
be fit by a power law or the fluxes extracted as described in Paper I. 



Eig. 1: The collection areas for standard and extended 
supercuts for the 1995/96 observing season are shown 
above. The resolution function is also plotted. 


3. METHOD 2: AN INTUITIVE APPROACH 

A different approach with emphasis on verifiability was developed at the University of Eeeds. There 
are two pieces to this approach, a method for selecting images likely to have been initiated by cosmic 
gamma rays and a method for determining the primary gamma-ray energy spectrum from the observed 
size spectra. Earlier descriptions of this approach are in (Hillas and West, 1991) and (West, 1994). 

The selection criterion is a “cluster” or “spherical” method in which a single parameter is used 
to characterize the gamma-ray-like nature of an image and correlations between image parameters are 
incorporated naturally. Simulated gamma rays produce images with four parameters (width, length, 
distance and alpha) that populate a four-dimensional space. Each real image can be tested to see if 
it is likely to be a gamma-ray image by the value of the Mahalonobis distance between it and the 
centroid of the cluster. (This is equivalent to scaling and rotation of the axes so that the window is 
spherical and fluctuations uncorrelated.) The selection window is defined by the expected position 
and dimensions of the gamma-ray parameter cluster for several broad ranges of image size. 




observations and the difference histogram is ascribed to gamma rays. A simulated size spectrum is 
then computed starting with a power-law primary energy spectrum for gamma rays. A weight is given 
to each simulated gamma ray and, by adjusting these weights, the spectrum can be varied so that its 
size spectrum matches that of the difference histogram. This method is simple, easy to implement and 
avoids the complexity of calculating collection areas and bias-free energy estimates. It is described in 
detail in Paper I. 

4. RESULTS 


The spectra obtained from the 1988/89 
and 1995/96 seasons using Method 1 and 
1988/89 season using Method 2 are in 
good agreement as shown in Fig. 2. As de¬ 
scribed in Paper I, we have tested the sen¬ 
sitivity of the results to uncertainties in the 
Monte Carlo simulations and have found 
that the results are relatively robust. The 
combined TeV fluxes from both seasons 
from Method 1 are well fit with the simple 
power law spectrum in which the differen¬ 
tial flux (J{E)) is given by: 

(3.3 ± .2 ± .7) 10^^ (_^)- 2 . 45 ±. 08±.05 

( 1 ) 

in units of m^^s^^TeV^^ where the first 
errors are statistical and the second are our 
estimate of systematic errors. 

However, the simple extrapolation of 
this fit to lower energies yields fluxes far in 
excess of those observed by EGRET as is 
clear from Fig. 3 which shows a quadratic 
fit of log(7) vs. log(£’) to our data and to an 
averaged point representing the EGRET flux (Nolan et ah, 1993) at 2 GeV. This fit may be written: 

J{E) = (3.25) jgY-i_ (2) 

5. COMMENTS ON INTERPRETATION 

Most of the early inverse Compton models for TeV gamma rays assumed a constant magnetic field in 
the principal source region where these are produced (e.g., Gould 1965 or Rieke and Weekes 1969) 
whereas more recent models (De lager and Harding, 1992 or Aharonian and Atoyan) incorporate 
hydrodynamic plasma/field flow making the calculations more complex and the results probably more 
realistic. Here, we try to stay close to the data and make the simpler assumption that the field is 
constant. The broad synchrotron emission band apparently extends up to 10^ eV and is boosted 
to higher energies via inverse Compton scattering. The scattering giving rise to TeV gamma rays 
occurs in the Klein-Nishina rather than in the Thomson scattering regime. This implies that the 
electrons giving rise to our detected gamma rays must have energies in the range of 2-10 TeV and the 
corresponding scattered photons would mostly have energies of about 0.005 to 0.3 eV. This conclusion 
is only very weakly model dependent (see Paper II). Since electrons with energies of a few TeV 
generate synchrotron radiation at about 0.4 keV in a field of about 25 nT (see next paragraph), the 



Eig. 2: The Crab spectrum in the range 0.3 to about 
8 TeV extracted using Methods 1 (open circles) and 2 
(solid circles) for the Whipple 1988/89 database and us¬ 
ing Method 1 (x’s) for the 1995/96 database are shown 
above. Also shown is a fit to the combined Method 1 re¬ 
sults (solid line) as well as an earlier spectrum (dashed 
line) taken from Vacanti et al. (1991). 



emitting TeV gamma rays. 

For assumed magnetic field values 
and the observed synchrotron flux, it is 
possible to deduce the spectrum of pri¬ 
mary electrons, presumably generated in 
the shock at the termination of the pul¬ 
sar wind (see, e.g., Coroniti and Kennel 
1985). From the ambient photon density 
and deduced electron spectrum, the TeV 
flux can be calculated, and results for B 
fields of 10, 20 and 40 nT are shown in 
Fig. 3. The shaded regions reflect un¬ 
certainties arising from ill-defined UV to 
soft-X-ray region of the synchrotron pho¬ 
ton spectrum. As can be seen from the fig¬ 
ure, the effective B field must lie between 
20 and 40 nT with 27 nT falling very near 
to our TeV data. Since even more ener¬ 
getic electrons only keep their energy for a 
short time, they should exist only near the 
pulsar wind shock. Hence, measurement 
of the TeV spectrum over a wider energy 
range may probe spatial variations in the 
nebular magnetic field (see Paper II). 
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Fig. 3: Whipple and EGRET observations of the Crab 
gamma-ray spectrum. The dotted line is a fit to the 
WHipple points and a single 2 GeVflux value. Full-line 
curves are predicted inverse Compton fluxes for 3 dijfer- 
enct assumed B fields. 



